Time-Symmetrized Kustaanheimo-Stiefel Regularization
نویسندگان
چکیده
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and weakly perturbed cases, numerical integration is necessary in situations where the perturbation is relatively strong. Kustaanheimo{Stiefel (KS) regularization is widely used to remove the singularity in the equations of motion, making it possible to integrate orbits having very high eccentricity. However, even with KS regularization, long-term integration is di cult, simply because the required accuracy is usually very high. We present a new time-integration algorithm which has no secular error in either the binding energy or the eccentricity, while allowing variable stepsize. The basic approach is to take a time-symmetric algorithm, then apply an implicit criterion for the stepsize to ensure strict time reversibility. We describe the algorithm in detail and present the results of numerical tests involving long-term integration of binaries and hierarchical triples. In all cases studied, we found no systematic error in either the energy or the angular momentum. We also found that its calculation cost does not become higher than those of existing algorithms. By contrast, the stabilization technique, which has been widely used in the eld of collisional stellar dynamics, conserves energy very well but does not conserve angular momentum. Subject headings: Numerical simulation, Regularization, Time-symmetric integration,Two-body problem
منابع مشابه
Efficient Orbit Integration by Kustaanheimo-stiefel Regularization Using Time Element
We have found that the introduction of the time element into the Kustaanheimo-Stiefel (K-S) regularization changes the rate of error growth from quadratic to linear with respect to the physical time, even when a scaling is not applied. It also mostly wipes out the periodic error components in the case of unperturbed orbits, independently of the application of the scaling. Since the increase in ...
متن کاملBasics of regularization theory
We consider the dynamics of three point masses, where we assume that the mass of the third body is so small that it does not affect the motion of the primaries. In the framework of the restricted three–body problem, we investigate the collisional trajectories, which correspond to a singularity of the equations of motion. We investigate the regularizing techniques known as Levi–Civita, Kustaanhe...
متن کاملReversible Adaptive Regularization: Perturbed Kepler Motion and Classical Atomic Trajectories
Reversible and adaptive integration methods based on Kustaanheimo-Stiefel regular-ization and modiied Sundman transformations are applied to simulate general perturbed Kepler motion and to compute classical trajectories of atomic systems (e.g. Rydberg atoms). The new family of reversible adaptive regularization (RAR) methods also conserve angular momentum and exhibit superior energy conservatio...
متن کاملReversible Adaptive Regularization I : Perturbed
Reversible and adaptive integration methods based on Kustaanheimo-Stiefel regularization and modiied Sundman transformations are applied to simulate general perturbed Kepler motion and to compute classical trajectories of atomic systems (e.g. Rydberg atoms). The new family of reversible adaptive regulariza-tion (RAR) methods also conserve angular momentum and exhibit superior energy conservatio...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 1996